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Abstract. We classified t>2 measurements according to their sensitivities w.r.t. 
to two planes, namely, reaction plane and participant plane. Likewise, in f2/e 
scaling, we showed that one needs to choose a e that is sensitive to the same 
jjj ■ plane as that u 2 is sensitive to. We presented our v^j t as a function of centrality 

and transverse momentum. We studied the ratio of v\lv\. We discussed the 
• applicable range for hydrodynamics, as well as implications to an incomplete 

^ ' thermalization. 
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1. Introduction 

(N 

qq , In ultra-relativistic heavy ion collisions, spectators pass through each other quickly. 

and the system begins its evolution with what is left behind - the overlap region 
of two nuclei. The pressure gradient convert the spacial anisotropy, quantified by 
eccentricity e, of the overlap region into anisotropy in momentum space, quantified 
by azimuthal anisotropy v 2 . t>2 is defined as the second Fourier coefficient in the 
description of particles distribution w.r.t. the reaction plane [Q], and it is largely 
determined by the collective motion in the in-plane direction, vi at Relativistic 
Heavy Ion Collider (RHIC) is reported to be large and, for the first time in heavy 
ion collisions, can be described by ideal hydrodynamics [[2]. Theoretical calculation 
shows that in order to explain the large V2 observed at RHIC, one has to assume 
that the shear viscosity is extremely small [|3]. That is one of the important reasons 
for which scientists think that a perfect liquid has been formed in relativistic heavy 
ion collisions [|4]. 

Since the announcement of the discovery of a perfect liquid, our understanding 
of the matter created at RHIC has continued to advance. Elliptic flow analyses 
have been extended to great details by all four experiments at RHIC. At the same 
time, due to different techniques used, it becomes an increasingly amount of work 
to understand/compare results across different experiments and, sometimes, even 
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different analyses within the same experiment. Therefore it becomes important to 
understand what is the relation between various vi measurements. For the physics 
side, in order to quantify how perfect the liquid is, it is necessary to re-examine 
the hydrodynamics limit. In this paper, we try to study the hydrodynamic behav- 
ior under a more general context, namely, the transport approach which recovers 
hydrodynamics when mean free path is extremely small if compared to the system 
size [[5] 



2. Choosing the right i> 2 and e pairs 

The ratio of t^/e reflects how well the initial anisotropy is converted into mo- 
mentum anisotropy [ [6] . This conversion process is directly affected by dynam- 
ics of the system, e.g. Equation Of State, thermalization etc. It is important 
to measure this quantity as accurate as possible. However, there exist many v 2 
measurements, for example, v 2 measured by event plane method (f2{EP}), by cu- 
mulants (v 2 {2} ^i^}) , by Lee- Yang zero method (-^{LYZ}), and by using event 
plane reconstructed with Shower Maximum Detectors at Zero Degree Calorimeters 
(^{ZDCSMD}), etc. Likewise there exist many e calculations. In this section, 
we will try to make connections between various v 2 and e methods, and make the 
justification for the right combination of them. 
Define 
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where er 2 = (x 2 ) — (x) 2 , a 2 — (j/ 2 ) — (y) 2 , and <r 2 y = (xy) — {y) (x), and the average is 
taken over the coordinates of the participants in a given event. With this definition, 
e x is the eccentricity of reaction plane (defined by the impact parameter), and e x 
is also called £pp. s y has a distribution centered at zero with finite width. The 
participant eccentricity measures the asymmetry in the participant plane (defined 

by the principle axis of the ellipsoid) , and is given by e par t = \j e2 x + £ y = £ pp 

When both e x and e y has a Gaussian distribution, to the first order this is 
true and it is supported by Glauber Monte Carlo simulations, then the probability 
density function for epp is 

dn e part ( £part(£Rp)\ ( £ lart + ( £ RpY \ , . s 

1 = ■ JL 2~Iq — 2 ex P — ^2 = B G{e part ; (e RP } , a e ), 

depart CT| V (T| J ^ 2(T| / 

(2) 

With this p.d.f., one can show that £ pa rt{4} = (srp) [Ej. Similarly, under the 
assumption that vi is proportional to the initial system eccentricity (this is not 
true over a broad centrality range, but for a fine centrality bin, it still ensures a 
good Gaussian for V2 and S2(= (sin2(</> — ^rp))), which is what is required for the 
derivation of BG formula shown in Eq.([5]) ), one can find that i>2{4} = (v2{RP}} [El 
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8 . That means that, v 2 {4} and £ pa rt{4} are measurements sensitive to the reaction 
plane, not the participant plane. Indeed, that explains the reason that STAR's 
w 2 {ZDCSMD} agrees with v 2 {4}, as shown in the left panel of Figffl v 2 {ZDCSMD} 
is measured with the first order event plane reconstructed by spectator neutrons 
thus is more sensitive to vi in the reaction plane (not participant plane). Because 
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Fig. 1. Left: u 2 measurements as a function of centrality. This plot is from 
Right: V2 of Au+Au collisions at y / Sj\T/v = 130GeV. This plot is made based on 
datapoints from [ [TO] 

the p.d.f. of the magnitude of the flow vector, the (/—distribution [HU], shares an 
almost identical formula as Eq. ([2]), the u 2 obtained from fitting the q— distribution 
should be equivalent to t> 2 {4}, as confirmed by experimental data in the right panel 
of FigUJ 

To summarize this section, we find that w 2 {4}, u 2 {ZDCSMD} and t> 2 {q — dist} 
are sensitive to anisotropy in the reaction plane, they should be scaled with the stan- 
dard eccentricity or 4-particle cumulant eccentricity. Other u 2 measurements that 
are based on two particle correlations, that includes v 2 {2}, w 2 {EP} and w 2 {scalarProduct} 
HU] , etc., should be scaled with participant eccentricity or the 2-particle cumulant 
eccentricity. 



3. Hints of Incomplete Thermalization 

The large data set from run IV Au+Au collisions at ^/snn — 200GeV allows us 
to extend w 2 {4} measurement to large pt and in fine centrality bins (see Fig. [2]). 
For the reason stated in the previous section, we scale i> 2 {4} by initial standard 
eccentricities. The left plot of Fig. [2] shows u 2 {4} scaled by the eccentricity from 
Color Glass Condensate (CGC) [[UJ, and the right plot shows i> 2 {4} scaled by the 
Monte Carlo Glauber eccentricity for wounded nucleons. As expected, for the CGC 
case, the magnitude of w 2 /e is lower if compared to the ratio in which a Glauber 
eccentricity is used. For both cases, we see the ratio rises from peripheral events to 
central events, indicating that stronger flow has been developed in central collisions. 
We also notice that the pt where u 2 reaches its maximum increases from peripheral 
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collisions to central collisions, which is consistent with the expectation that the 
applicable range for hydrodynamics extends to large pt in central collisions. Note 
that vijt shows sign of saturation in central collisions for the CGC case but not 
much for the Glauber case. This is explained by [ [11] as the following: in very 
peripheral collisions, due to little asymmetry in the saturation scales, the CGC 
eccentricity approaches the same value as that in the Glauber model, but in central 
collisions, CGC predicts a larger eccentricity than the Glauber model when there 
is a large asymmetry in the local saturation scales of the collisions partners, along 
a path in impact-parameter direction away from the origin. Note that i>2{2} is not 
suitable for this study because it is more susceptible to nonflow at large pt- 
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Fig. 2. V2 scaled by initial CGC eccentricity (left) and Glauber eccentricity (right) 
as a function of pt- This plot is from [H"2]. 

To understand how well hydrodynamics describes STAR's V2, we investigated 
the behavior of V2 under the contex of the transport model, which will be reduced 
to hydrodynamics when the mean free path is much smaller than the system size [ 
5 . In such approach, the dependence of V2/e on particle's density in the transverse 
plane (1/SdN/dy) can be described by: 

(3) 



f'2 

e 



hydro 1 + K/K a 



hydro ^ j_ j ^c^^dN] J_ 



c 4S dy I K, 



Where K is Knudsen number defined by the mean free path divided by the system 
size (sometimes it is more convenient to use K~ 1 which means number of collisions 
a particle encounters before it escapes), and Kq is a constant can be determined 
through transport calculations. In our study, we take K = 0.7 following the sugges- 
tion of Ollitrault [[11]. A factor of 4 in front of S is to take into account the different 
definition of S between STAR (S = tt^/{x 2 ) (y 2 )) and [E]. ((S = ^^(x 2 ) (y 2 )). 
\re'\ hydro anc ^ a are ^ ree P arame ters that have to be determined from fitting the 
data. In this approach hydrodynamic limit of i»2/e can be never reached, but can 
be only asymptotic approached. 

Fig. [3] shows «2{4} scaled by CGC initial eccentricity and Glauber initial ec- 
centricity for wounded nucleons, as a function of particle density in the transverse 
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f*ig. 3. W2{4} scaled by initial eccentric- 
ties, as a function of centrality. Fitted 
lydrodynamic limits for v^jt are indi- 
Fined'Lim'itwitKcGc'e oated by horizontal lines. CGC eccen- 
;ricity and overlap area S are from [ lllj . 
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plane, with fits to Eq. [3l The fitted hydrodynamic limits for the ratio of vi /e are 
0.23 and 0.36, for CGC case and Glauber case, respectively. For the same reason 
mentioned above, we see that the curve shows a hint of saturation for the CGC case, 
but not for the Glauber case, due to the relatively larger e for CGC case in central 
collisions. It is interesting to see that, for central Au+Au collisions, the ratio of 
i>2/e is about 20 — 30% away from hydrodynamic limits. It means that, there is 
still significant room for flow to grow before the system saturates at hydrodynamic 
limits. 

From the simple observation that both V2 and V4 are proportional to A -1 
for small A -1 , one expects that Vi/v^ decreases with A -1 , reaching a minimum 
when the hydrodynamical regime is reached. For this reason the ratio of v^jv\ has 
been argued as a probe to test the degree of thermalization. Fig. [4] shows STAR's 
measurement of v^jv\ as a function of transverse momentum. The major systematic 
uncertainty in this measurement comes from vi [[12]. In this analysis, the induced 
systematic error from V2 uncertainty is estimated by studying the difference between 
i>2{4} and V2 measured with event plane constructed by tracks from STAR's Forward 
Time Projection Chamber (FTPC). This advanced study reduces the previously- 
reported systematic error at QM06 conference by 40% relatively. The dashed lines 
are ratio come out of calculations by solving Boltzmann equations with Monte 
Carlo simulation, with different Knudsen number A. When the Knudsen number is 
small, it recovers the hydrodynamic limit as indicated the solid line. The plot shows 
that the system exhibits, again, significant deviation from ideal hydrodynamic limit 
(A" << 1), and the data is consistent with a incomplete thermalized system with 
A > 0.5. 



4. Conclusion 

To summarize, we found that i>2{4}, ^{ZDCSMD} and i>2{q— dist} are all sensi- 
tive to azimuthal correlation w.r.t the reaction plane, not the participant plane, thus 
they should be scaled by eccentricities that are sensitive to reaction plane too. That 
includes standard eccentricity and 4-particle cumulant eccentricity. For vi methods 
that are based on two particle correlations, they are sensitive to the azimuthal cor- 
relation w.r.t the participant plane, and they need to be scaled by the corresponding 
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Fig. 4. W4{EP 2 }/w|{4} as a function of 
p t - This plot is from [[12] with an ad- 
vanced estimation on systematics repre- 
AuAu 200 GeV 20 - 30%scnted by half box-brackets. 
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eccentricities that are sensitive to the participant plane. That includes participant 
eccentricity or 2-particle cumulant eccentricity. We found that from peripheral to 
central Au+Au collisions flow increases, and the applicable range for hydrodynam- 
ics extends to larger p t . However, V2/e and v^jv\ shows significant deviation from 
ideal hydrodynamic limit, when that limit is extracted from fitting the data itself 
with a Boltzmann equation motivated formula. Our study shows that although in 
general hydrodynamic does a good job in terms of describing w 2 at RHIC, there are 
features that are not consistent with a complete thcrmalization and they cannot be 
easily dismissed. 
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